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ABSTRACT 

Medvedev & Narayan have described a hot accretion flow onto a spinning neutron star 
in which the gas viscously brakes the spin of the star. Their self-similar solution has 
the surprising property that the density, temperature and angular velocity of the gas 
at any radius are completely independent of the outer boundary conditions. Hence, the 
solution cannot be matched to a general external medium. We resolve this paradoxical 
situation by showing that there is a second self-similar solution which bridges the gap 
between the original solution and the external medium. This new solution has an extra 
degree of freedom which permits it to match general outer boundary conditions. We 
confirm the main features of the analytical results with a full numerical solution. 



1 INTRODUCTION 



Accretion flows around compact objects frequently radi- 
ate significant levels of hard X-rays. This has provided 
strong motivation fo r the study of hot accretion flows . 
IZeldovich fc Shakural il969f) and lAlme fc Wilson] dl973l) 
considered spherically free-falling plasma impinging on the 
surface of a neutron star (NS). They calculated the pene- 
tration depth of the falling protons and made preliminary 
estimates of the radiated spectrum. Their i deas were fol- 
lowed up by a number of later authors (e. g., iTurolla et alJ 
ll994l:IZampieri et alll995HZane et alll998l). w ho com puted 
detailedspectran^^ (l200lT) modi- 

fied the model of lZeldovich fc Shakural B1969I) by considering 
a rotating advection-dominated accretion flow around a NS. 
Their model represents an improvement on the earlier work 
since it includes angular momentum and viscosity. 

A fluid approach to hot ac c retion onto a NS was pi- 
oneered by IShapiro fc Salpeterl <ll975l) . who worked out 
the structure of the standing shock in a spherical flow, 
and computed the two-temperature structure and the re- 
sulting radiation spectrum of the post-shock gas. The 
equivalent proble m for an a c cretin g white dwarf was 
analyzed bv [ K vlafis fc LambI (|W82!). In related work, 
IChakrabarti fc Sahul J1997F described the hydrodynamics of 
spherical accretion onto black holes and NSs, but without 
including radiation processes. 

The above studies involve flows in which the ac- 
creting matter crashes on the surface of the star, form- 
ing a discontinuity or a s hock of some kind. Recently, 
IMedvedev fc Naravar] fcOOll hereafter MNOl) discovered a 
rotating solution of the viscous fluid equations that corre- 
sponds to hot quasi-spherical accretion onto a spinning NS. 
Their solutio n is closely related to the two-temperat ure flow 
described bv IShapiro. Lightman. fc EardlevI il976T> . A fea- 
ture of the MNOl solution is that the gas moves subsoni- 



cally in the radial direction and merges with the accreting 
star without a shock. The flow essentially "settles" onto the 
rotating star; the solution may thus be referred to as a "hot 
settling flow" (hot because the gas is at the virial tempera- 
ture and has a quasi-spherical morphology). MNOl showed 
that the accreting gas removes angular momentum from the 
central star and that this braking action dominates the en- 
ergy equation of the accreting gas. The flow could thus also 
be called a "hot brake." The hot settling flow should not 
be confused with th e boundary layer which for ms close the 
stellar surface (e.g., iNaravan fc Pophamlll99'3) . where the 
gas density is high and steep spatial gradients are present. 
The hot settling flow forms outside the boundary layer and 
extends radially to a large dista nce, typically thousan ds of 
stellar radii or more (see MNO l) . Ilkhsanovl fcOOll l2003ft . fol- 
lowing up on earlier work bv iDavie^^^Prjngla fl98l"l) . has 
recently described a subsonic hot accretion flow around a 
magnetized neutron star in the propeller state. The relation 
between his solution and our hot settling flow is discussed 
in §4. 

The relevance of the hot settling flow to real sys- 
tems is presently unclear, though several accreting white 
dwarf and black hole systems have bee n suggested as candi- 
dates for such a rotation-powered flow jMedvedev fc Menoul 
12001 IMedvedev fc Murravl 120021) . While MNOl described 
the properties of the self-similar region of the flow, they did 
not discuss how to match the solution to realistic boundary 
conditions. The matching to an external medium at large 
radii is particularly problematic, since the self-similar set- 
tling flow solution has the remarkable property that the 
density, temperature and angular velocity of the gas at given 
radius depend only on the dimensionless spin of the central 
star, and are completely independent of the properties of 
the external medium. This leads to an apparently serious 
problem. If one tries to match the solution to the external 
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medium by selecting the radius at which the pressure of the 
solution matches that of the medium, then neither the den- 
sity nor the temperature will agree. Does this mean that 
the solution is physically inconsistent? We argue otherwise 
in this paper. 

We show that there exists a second self-similar solution 
at radii outside the MN01 hot settling flow, which acts as a 
bridge between the MN01 solution and the external medium. 
This bridging solution has an extra degree of freedom which 
allows it to match a general density and temperature in the 
external medium. Apart from the additional degree of free- 
dom, the solution retains many of the features of the MN01 
settling flow: (i) it is pressure supported, (ii) it resembles a 
static atmosphere at low mass accretion rates, and (iii) the 
angular momentum flux associated with the viscous braking 
of the central star dominates most of the physics. 

The paper is organized as follows. In J^we derive ana- 
lytical expressions for various self-similar solutions. In §^]we 
present full numerical solutions and compare them with the 
analytical solutions. Finally, in §4 we conclude with a brief 
discussion. 



2 ANALYTICAL SELF-SIMILAR SOLUTIONS 

We consider gas accreting viscously onto a compact spin- 
ning object. The central object has a radius R,, a mass 
M* = toMq, and an angular velocity = sQk{R*), where 
Qk{R) = (GM/R 3 ) 1/2 is the Keplerian angular velocity 
at radius R. We measure the accretion rate in Eddington 
units, rh = M/M-^dd, and the radius in Schwarzchild units, 
r = R/Rg, where A-fEdd = 1.39 x 10 18 m g s _1 (corresponding 
to a radiative efficiency of 10%) and R g = 2GM/(? . 

We assume that the flow is hot and quasi-spherical, 
which generally requires a low mass accretion rate (see 
Narayan et al. 1997). The accreting gas has nearly the virial 
temperature, i.e., c 2 ~ GM/R ~ (QkR) 2 , and the local 
vertical scale height H = c 3 /Q,k is comparable to the local 
radius R. We may then use the height integrated hydrody- 
namic equations for a steady, rotating, axisymmetric flow, 
and for simplicity we may set H — R (see MN01). (Note 
that, even when H m R, the stability properties of the flow 
may depend on whether we use H or R in the equations; 
this is briefly discussed in §4). In the following, we closely 
follow the analysis of MN01, with a few changes. 

For simplicity, we set rh = and omit the continu- 
ity equation. Thus, the gas configuration corresponds to a 
radially static "atmosphere." The motivation for this ap- 
proximation follows from the observation that the density 
p, temperature T and the angular velocity Q of the gas in 
the MN01 self-similar solution are completely independent 
of rh. Only the radial velocity v depends on rh, and it is 
given trivially by the spherical continuity equation 



A I 



4nR 2 p' 



(1) 



Because of this, we do not lose any generality by setting 
rh = v — in the analysis; we may always introduce a finite 
rh and finite v after the fact. 

We assume that the flow is highly sub-Keplerian, 
Q(R) <C Q,k(R), so that the centrifugal support is negligible 



compared to the pressure support. The radial momentum 
equation then takes the following simple form, 

GM _ 1 djpcp 

and written the 



(2) 



p dR ' 

where we have used the fact that v 
pressure as p = pc 2 where c s is the isothermal sound speed. 
MN01 present a more complete analysis in which they do 
not assume that the rotation is slow. They then obtain an 
extra factor of (1 — s 2 ) in their equation, which propagates 
through to all the results. Since we ignore the factor, our 
analysis corresponds to the case of a slowly-spinning star: 
s 2 < 1. This approximation is made only to simplify the 
analysis, and all the results may be generalized for arbitrary 
s. 

From the analysis in MN01, we know that the accreting 
gas in our problem acts as a brake on the central spinning 
star and transports angular momentum outward through the 
action of viscosity. We therefore write the angular momen- 
tum conservation equation for the gas as follows, 

, 4 dfi 



J = A-KvpR" 



dR 



constant, 



(3) 



where j is the outward angular momentum flux, and v is the 
kinematic coefficient of viscosity. This equation is exactly 
valid in steady state if rh = 0. When rh is non-zero, there is 
an additional term, MflR 2 , due to the flux of angular mo- 
mentum carried in by the accreting gas. The key feature of 
the MN01 hot brake solution is that the latter flux is negli- 
gible compared to the outward flux from the star. Equation 
(j^J is, therefore, valid even when rh 7^ 0, so long as rh is 
small enough for the term MQR 2 to be negligible. 

We employ the usual a prescription for the kinematic 
coefficient of viscosity, which we write as 



v = ac a H « ac s R- 



(4) 



Often, in accretion problems, one makes use of the rela- 
tion H = Cs/Qk and writes v = ac 2 /fijc. This prescription 
is equivalent to equation (4) in the regime of the MN01 
hot settling flow. However, in the outer regions of the flow, 
where the two new solutions described in §§2.2,2.3 apear, H 
is much less than c s /Q,k, and v — ac 2 /Q,K is not a good ap- 
proximation. Equation (4) is a superior prescription and is 
physically better motivated over a wide range of conditions 
(so long as the flow is quasi-spherical). 

For simplicity, we assume that the gas is one- 
temperature; it is straightforward to generalize to the 
two-temperature regime, as done in MN01, but the one- 
temperature analysis suffices for the present paper. Hence 
do not need to treat electrons and protons separately. 

Viscous braking heats the accreting gas, and we assume 
a steady state in which this heating is balanced by cooling. 
The energy conservation equation thus becomes 

q + = q~, (5) 

where the viscous heating rate per unit volume q + and the 
radiative cooling per unit volume via bremsstrahlung q~ are 
given by (MN01) 
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where a/ is the fine structure constant, <tt is the Thomson 
cross-section, and m e and m v are the mass of the electron 
and the proton. Equation ® is exactly valid for a steady 
state flow with m = v — 0. When v 7^ 0, the energy equation 
has another term corresponding to the advection of energy. 
In advection-dominated accretion flows, for instance, this 
term dominates over the cooling term q~ jNaravan et alJ 
119971) . In the present case, however, we consider a situation 
in which the advection term is negligible (which corresponds 
to low rh). 

Finally, we assume that the spinning star is immersed in 
a uniform external medium with a density p ex t, temperature 
T cx t and pressure p ex t- We seek an accretion flow solution 
that extends from the spinning star on the inside to the 
external medium on the outside. As we show below, the 
solution consists of two distinct self-similar regimes, plus a 
third asymptotic regime inside the external medium. 

2.1 Inner self-similar solution 

We first consider the inner regions of the flow, where the 
pressure p p ox t- This is the regime of the MN01 solution, 
where the variables have the following radial dependences: 



P = pir 



T = Tir" 



(6) 



The subscript "1" in the coefficients is to indicate that this 
is the first solution, to distinguish it from the second and 
third solutions described below. By substituting the above 
solution in equations © and JSJ, we see that it satis- 
fies the basic conservation laws. We may also solve for the 
numerical constants: 



Pi = 

fcTi = 
Qi = 



as 2 9 
"r7 W> 



1/2 
m e \ 3 

C <Vff,NR, 

m„ ■ 
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s c 



= S Q K (Rg 



V2R g 

We note that if rh 7^ then the flow has a small constant 
radial velocity: 

v oc r , (7) 

as follows from equation |T). 

The angular momentum flux in the solution is given by 



J 



o5/2 
2 3 r>2 o 



-a s R, 



9 2 5/2 
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Qff.NR 



(8) 



By assumption, this flux is much greater than the angular 
momentum flux due to accretion, which sets an upper limit 
on the mass accretion rate for the solution to be valid (see 
MN01). The pressure is given by 
-3 



V 



2 2—3 

pc s = p!C sl r 



Pir, (9) 

where <? sl — 2kT\/m p . 

The above self-similar solution describes the flow at 
radii r <C (pi/p cx t) 1//3 , where the pressure p 3> p cxt . As 
mentioned in §1, the solution has the remarkable property 
that all the quantities are uniquely determined by a single 
parameter s — the dimensionless spin of the central ob- 
ject — specified on the inner boundary. The fact that the 
solution does not depend on the outer boundary condition 



in any way means that there is no simple way to match it 
to the external medium. Clearly, there has to be a second 
solution to bridge the gap between this solution and the ex- 
ternal medium. We derive the bridging solution in the next 
subsection. 



2.2 A second self-similar solution 

We consider next the gas that lies just outside the region 
of validity of the first self-similar solution described above. 
In this zone, the pressure is expected to be approximately 
equal to the external pressure p ex t: 



pc s = pext = constant. 



(10) 



This condition replaces the hydrostatic equilibrium equation 
while equations ||3J and (J^J continue to be valid. In this 
region, we find that there is a second self-similar solution of 
the form 



P = P2T 



-7/2 



T = T 2 r 



7/2 



(ii) 



where the label "2" refers to the fact that this is our second 
solution. 

To match the second and first solutions, we require 
that the fluxes of angular momentum in the two solutions 
must be equal; this yields the constraint (3/2)pifiiT 1 1//2 = 
(9/4)p 2 n 2 T 2 1/2 . Making use of this and the other equations, 
we solve for the numerical coefficients in equation 1111 : 



P2 = 



kT 2 = 



a 3 /V 3 5/2 / m E 



pI^rT 29/4 \«h> 



3/4 



Vff.NR 
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Q 3/2„3 



v 1/4 R 5/4 
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The pressure in this solution is constant and equal to the 
external pressure, p GX t, and the angular momentum flux is 
also constant and is equal to J in equation 0. If the flow 
has a small but nonzero accretion rate, m/0, then its radial 
velocity varies as [see eq. ijljl] 



3/2 



(12) 



Whereas the original MN01 self-similar solution has a 
unique profile for a given choice of s, we see that the sec- 
ond solution derived here has an extra degree of freedom, 
namely the external pressure p ox t . This extra degree of free- 
dom solves the problem discussed in §1. Thus, the full so- 
lution consists of two zones: an inner zone described by the 
first (MN01) solution © and an outer zone described by 
the second solution Qjl . The radius r ma t c h at which the 
two solutions match is obtained by equating the pressures: 



W 8 3 1 / 3 fm e Y /G 



f"match 



S/3 n l/3 



(13) 



2^6 y mp/ 

The second solution matches the external medium at 
the radius r cx t at which its temperature matches that of the 
medium. This gives 

3/14 



3 5/7 



r ex t 



8/7 



•,3/7 



(14) 
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If we wish we could also write this in terms of the external 
density by making the substitution fcT oxt = m p p ox t/2p e xt- 

2.3 Asymptotic solution in the external medium 

For completeness, we present here the solution inside the ex- 
ternal medium. By assumption, the external medium has a 
uniform temperature and density, and a uniform rate of cool- 
ing. To maintain equilibrium, there has to be some constant 
source of heat that exactly compensates for the cooling. We 
assume that such a source of heat exists (e.g., cosmic rays). 
The rotation Q is non-zero, but it decays rapidly outward. 
The small amount of rotation helps to transport the angular 
momentum flux from the star out into the external medium. 
Solving the angular momentum conservation law (3), we ob- 
tain the following solution 

p = p cx t, T = T cx t, = f2 3 r~ 4 , (15) 
where 

" s = a * 3 $^iwr-p^ (^ C xt)- 1/2 , (i6) 

2'tt KgQs^B. 

and pcxt = 2fcT xtPext/?Tip. For m / 0, the velocity scales as 
r~ 2 . 



3 NUMERICAL RESULTS 

The three self-similar solutions written above are special so- 
lutions of the basic differential equations (2), (3) and (5), 
which are valid under specific conditions. To check the va- 
lidity of these analytical solutions, we have computed numer- 
ical solutions of the basic differential equations. We use the 
same code as in MN01, with two changes. First, we switched 
to the viscosity prescription given in equation @, rather 
than the prescription v = ac? a /Q,K used in MN01. Second, in 
addition to viscous heating, we included a constant heating 
rate which we adjusted so as to balance the radiative cooling 
in the homogeneous external medium (see §2.3). The code 
uses a relaxation method to solve the one-dimensional hy- 
drodynamic equations with specified inner and outer bound- 
ary conditions. Although it employs the full equations of a 
two-temperature plasma, the results are essentially equiva- 
lent to those of a one-temperature plasma in the region of 
interest for this paper, namely the region at large radius 
where the flow matches onto the external medium 

In the calculations, the flow was taken to extend from 
an inner radius i?i n = 3 R g to 7? ou t = 10 7 R g . The mass ac- 
cretion rate was taken to be low, m = 2 x 10 -5 , in order that 
the flow should correspond to the regime of the hot settling 
flow solution. We took the viscosity parameter to be a = 0.1 
and set the spin of the star to be s — 0.3 (i.e., 30% of the 
Keplerian rotation at the stellar surface). We took the other 
inner boundary conditions to be the same as in MN01. At the 
outer boundary, we specified the temperature and density of 
the external medium. Figure^shows four solutions. The ex- 
ternal temperature is kept fixed at T(i? C xt) = 10 8 K in all the 
solutions, but the external density varies by a decade and a 
half: p(Rcxt) = 2.5 xlO 9 , 8.1x10 s , 2.5x10 s , 8.1xlO r cm~ 3 . 
We have also done other calculations in which we kept p ex t 
fixed and varied T ex t • These give very similar results. 

Fig. shows that, right next to the star, there is a 



boundary layer, where the density rises sharply as one goes 
into the star and the temperature drops suddenly. We do 
not analyze this region. Once we are outside the bound- 
ary layer, the gas behaves very much according to the an- 
alytical solutions discussed in §2. Starting just outside the 
boundary layer and extending over a wide range of radius, 
the numerical solution exhibits a self-similar behavior with 
power-law dependences of the density, temperature and an- 
gular velocity. This region corresponds to the self-similar 
solution of MN01. There are, in fact, two zones, an inner 
two-temperature zone, and an outer one-temperature zone 
(MN01). The latter corresponds to solution 1 (eq. dis- 
cussed in §2.1. The most notable feature of this region is 
that the density, temperature and angular velocity of the 
numerical solutions are completely independent of the outer 
temperature and density, as predicted by the analytical so- 
lution. The slopes of the numerical curves also agree well 
with the analytical scalings. 

At a radius i? matc h ~ 5 x 10 4 ...2 x 10 5 R g 9depend- 
ing on the outer pressure, see eq. 1131 . solution 1 merges 
with solution 2 (eq. II II described in §2.2. Here, the solu- 
tion does depend on the outer boundary conditions, and it 
scales roughly according to the slopes derived analytically. 
At even larger radii R > f? cx t ~ 3 x 10 5 ...2 x 10 6 R g (see eq. 
I14|l . the flow matches onto the ambient external medium. In 
this region we have solution 3 (eq. H5H described in §2.3. As 
expected, out here only the angular velocity and the radial 
velocity vary with radius. Both have the scalings predicted 
for solution 3. 



4 DISCUSSION 

In this paper, we have removed one piece of mystery sur- 
rounding the self-similar "hot settling flow" or "hot brake" 
solution discovered by MN01. Specifically, we have shown 
that the remarkable insensitivity of the MN01 solution to 
external boundary conditions is a consequence of the fact 
that the solution is insulated from the outer boundary by 
the presence of a second solution, which bridges the gap 
between the first solution and the external medium. We de- 
rived the form of the second solution analytically in §2.2 
and showed via numerical computations (§3, Fig. 1) that 
the two solutions together are able to match a wide range of 
outer boundary conditions. This solves one of the mysteries 
associated with the hot settling flow solution. 

There are, however, two other problems that still need 
to be addressed. First, the solution we have derived treats 
the mass accretion rate m as a free parameter. (Indeed, the 
analytical solutions were obtained for the limit m — > 0, i.e. 
for a hot atmosphere.) What determines rh? It is certainly 
not the outer boundary, since we have obtained the com- 
plete outer solution. The accretion rate must therefore be 
determined by an inner boundary condition. This is not un- 
expected. In the case of spherical accretion, one recalls that, 
while the accretion rate for the transonic solution is deter- 
mined by the outer boundary conditions, the accretion rate 
for the subsonic settling solution is determined by an in- 
ner boundary condition. In that problem, a whole family of 
settling solutions exists. Each member of the family has a 
different mass accretion rate, and it is the manner in which 
the gas cools and condenses on the accreting star that deter- 
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mines which particular solution, i.e., which rh, is selected. 
We expect the same situation to apply to our problem. Un- 
fortunately, this means that in order to estimate rh we have 
to solve the full coupled hydrodynamic and radiative trans- 
fer equations in the boundary layer region next to the neu- 
tron star, with proper boundary conditions. We have not yet 
succeeded in this difficult exercise. 

The second problem that needs to be addressed is the 
stability of the hot settling flow solution. Because the solu- 
tion is hot, optically thin and satisfies the energy balance 
condition it is similar in many respects to the hot solu- 
tion discovered by Shapiro et al. (1 976). The lat ter is known 
to be thermally very unstable (e.g.,|Piranl ll978f) . and so one 
wonders whether the hot settling flow might also be un- 
stable. The thermal stability depends, in particular, on the 
viscosity prescription chosen, as illustrated below (we thank 
the referee for this simple argument). With the prescription 
v = ox? s /Q.k the flow is unstable because q + oc T while 
q~ oc yT. Hence, a local increase of temperature results in 
an increase of the net heating rate q + — q~ which leads to a 
further temperature increase. However, for the viscosity pre- 
scription given in equation (3) , we have q + oc vf oc q~ , so 
that the flow is marginally stable. A more detailed analysis 
of the problem is beyond the scope of the present paper. 

We finally comment on the relationship of the hot set- 
tling flow discus sed in this paper to th e subsonic propeller 
flow described by I khsanovl il2()()lL|2(lll. i : . Both flows describe 
the braking action of hot gas on a spinning neutron star. 
The propeller flow has been discussed in connection with 
a strongly magnetized neutron star, while the hot settling 
flow was developed to model accretion onto an unmagne- 
tized neutron star, but this is not a large distinction. The 
main difference between the two solutions is in the treat- 
ment of the energy equation. In fkhsanov's subsonic pro- 
peller flow, the heating rate of the accreting gas through 
viscous dissipation is much larger than the radiative cooling 
rate. The gas becomes convective and isentropic, with den- 
sity falling as r _3//2 . The solution is in some sense related 
to an advection-dominated accretion flow (Narayan et al. 
1997) or a convection-dominated accretion flow (Narayan, 
Igumenshchev & Abramowicz 2000; Quataert & Gruzinov 
2000). In contrast, the hot settling flow, as well as the two 
other solutions described in this paper, satisfy detailed en- 
ergy balance at each radius. The viscous heating at each 
point is exactly balanced by local radiative cooling via opti- 
cally thin bremsstrahlung, as indicated in equation (5) of the 
present paper. Furthermore, the density falls off as r" 2 and 
r~ 7 / 2 for solutions 1 and 2 rather than r~ 3 ^ 2 , the entropy of 
the gas increases outward, and the gas is convectively stable 
(MN01). 

The authors are grateful to the referee for helpful sug- 
gestions, and RN thanks Lars Hernquist for useful discus- 
sions. This work was supported in part by NASA grant 
NAG5-10780. 
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Figure 1. Profiles of density (top left panel), temperature (top right panel, the electron temperature is the lower curve on the left and 
the proton temperature is the higher curve), angular velocity (bottom left panel), and radial velocity (bottom right panel), for four 
numerical solutions of the full height-integrated differential equations. The four solutions correspond to different values of the density 
of the external medium: p C xt — 2.5 x 10 9 , 8.1 x 10 s , 2.5 x 10 8 , 8.1 x 10 7 cm . The first and fourth solutions are labeled 1 and 4, 
respectively. The temperature of the external medium and the accretion rate are kept fixed in all the solutions: Tp ext = T e oxt = 10 8 K, 
m = 2 X 10~ 5 . The analytical slopes of the three self-similar solutions described in S H2.1ll2~^l are shown for comparison 



